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1 | INTRODUCTION 


Bovine coronavirus (BCoV) is a causative agent of respiratory and 
enteric diseases in calves and adult cattle. It belongs to the virus 
species Betacoronavirus 1 of the lineage A of the genus Betacoron- 
avirus (Bidokhti et al., 2013). 

BCoV has relatively recent common ancestors with human coro- 
OC43 (HCoV-OC43) 
encephalomyelitis virus (PHEV) (Vijgen et al., 2006). Due to their 


genetic closeness, the interspecies transmission potential of BCoV 


navirus and porcine haemagglutinating 


should not be ignored. In fact, human enteric coronavirus 4408 
strain (HECV-4408) was first isolated from a child’s diarrhoeic stool 
and was genetically and antigenically more closely related to BCoV 
than HCoV-OC43 (Zhang, Herbst, Kousoulas, & Storz, 1994). 


*These authors contributed equally to this work 


Bovine coronavirus (BCoV) is a causative agent of respiratory and enteric diseases 
in cattle and calves. BCoV infection was also evident in captive wild ruminants. 
Recently, water deer are recognized as the most common wildlife to approach farm- 
houses and livestock barns in Korea. Therefore, we investigated 77 nasal swab sam- 
ples from non-captive wild water deer (Hydropotes inermis) between November 
2016 and September 2017 and identified three samples positive for coronavirus, 
indicating potential for respiratory shedding. The full genomic sequences of the 
water deer coronavirus were closely related to BCoV (>98%). Therefore, effective 
biosecurity system in bovine farms would be necessary to prevent contact between 


farm ruminants and free-ranging wild water deer. 


bovine coronavirus, carrier, nasal swab, water deer 


In the late 1970s, bovine-like CoVs were first detected in other 
ruminants by electron microscope (EM) and they have since been 
detected in the faecal samples of wild ruminants such as sambar 
deer (Cervus unicolor), waterbuck (Kobus ellipsiprymnus), sable ante- 
lope (Hippotragus niger), white-tailed deer (Odocoileus virginianus), 
wisent (Bison bonasus), Himalayan tahr (Hemitragus jemlahicus), sita- 
tunga (Tragelaphus spekii), nyala (Tragelaphus angasii) and giraffe 
(Giraffa camelopardalis) (Alekseev et al., 2008; Chung, Kim, Bae, Lee, 
& Oem, 2011; Hasoksuz et al., 2007). As no specific genetic markers 
were identified to discriminate the ruminant bovine-like CoVs from 
BCoVs (Alekseev et al., 2008), interspecies transmission of BCoVs 
may be common among the ruminants. 

However, the majority of cases were found in the captive wild 
ruminants rather than non-captive ruminants and most of them were 
detected in the faecal samples. In Korea, no evidence of BCoV infec- 


tion in water deer (Hydropotes inermis) was reported between 2010 
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TABLE 1 Information on nasal swab samples 





Month of Collection 


November December January February March 
2016 2016 2017 2017 
No. of nasal 1 8 5 2 5 
swabs 
No. of CoV O 1 1 1 O 
positives? 


2017r 


KIM ET AL. 
April May June July August September 
2017 2017 2017 2017 2017 2017 Total 
10 11 16 7 7 5 77 
O O O O O O 3 


“Positives in RT-PCR based on consensus primers targeting RNA-dependent RNA polymerase region of coronaviruses (Poon et al., 2005). 


and 2012 (Kim et al., 2014). Therefore, we investigated coron- 
aviruses in water deer in Korea, using 77 nasal swabs collected 
between November 2016 and September 2017. 


2 | MATERIALS AND METHODS 


2.1 | Samples 


In this study, 77 nasal swab samples from water deer (Hydropotes in- 
ermis) were obtained from Chungnam Wild Animal Rescue Center, 
between November 2016 and September 2017. The water deer in 
this study were rescued wildlife, and the nasal swabs were immedi- 
ately collected as soon as they were delivered by wildlife veterinari- 
ans. The nasal swab samples were placed into universal viral 
transport medium (BD Diagnostics, USA) and transported to the lab- 


oratory for analysis. 


2.2 | Viral detection 


RNA was extracted using TRIzol LS (Invitrogen, USA), and consensus 
primer-based reverse transcriptase polymerase chain reaction (RT- 
PCR), based on the RNA-dependent RNA polymerase (RDRP) gene, 
was performed to detect coronaviruses (Poon et al., 2005). Among 
77 samples, three samples appeared positive from RT-PCR and the 
amplicons were sequenced using target-specific forward and reverse 
primers synthesized by Cosmo Genetech Co., Ltd. The nucleotide 
sequences were further analysed against other coronavirus 
sequences from GenBank using BioEdit (http://www.mbio.ncsu.edu/ 
BioEdit/bioedit.html) and MEGA version 6 (Tamura, Stecher, 
Peterson, Filipski, & Kumar, 2013). 


2.3 | Genome sequencing 


The genome of water deer coronavirus was further sequenced by 


RT-PCR using primers that were designed based on the genome 


of a related BCoV, strain Quebec (AF220295). Subsequently, 3' 
and 5' end sequences were obtained by rapid amplification of 
cDNA end (RACE) sequencing based on PCR using the adapter- 
oligo dTVN primer from Cosmo Genetech Co., Ltd. The obtained 
partial viral genome was assembled and manually annotated by 
comparisons with the open reading frames (ORFs) of other CoVs 
in GenBank database and further schematized with Geneious Ro 
(Biomatters Ltd., Auckland, New Zealand). To compare the com- 
plete genome of the water deer coronavirus and other related 
CoVs, the orthologous average nucleotide identity (OrthoANl) 
value was preliminarily analysed using ANI calculator (http:// 
www.ezbiocloud.net/tools/ani), and the nucleotide and deduced 
amino acid sequences of the viral ORFs were compared by multi- 
ple sequence alignment using CLUSTALW software (www.ebi.ac. 


uk/clustalw), respectively. 


3 | RESULTS AND DISCUSSION 


The three positive samples (W17-12, W17-15 and W17-18) were 
obtained in December 2016, January 2017, and February 2017, 
respectively, which covered the winter season in Korea (Table 1). 
Phylogenetic analysis, based on 256 bases of the partial RDRP gene, 
showed they were closely related to BCoV and HECV-4408 (Fig- 
ure 1a). The nucleotide identities of the viruses with BCoV and 
HECV-4408 were 99.2% and 100%, respectively. 

The fully sequenced water deer coronavirus W17-18 genome 
(GenBank Accession No. MG518518) was 31,034 bp (37.1% G+C 
content); the overall gene content and arrangement in the strain 
were almost identical to those from its related BCoVs from captive 
wild ruminants in the GenBank database. Moreover, the OrthoANl 
value analysis suggested that the water deer coronavirus W17-18 
genome was very similar (>98%) to those BCoVs from bovine and 
wild animals and closely related to HCoV-OC43 strains (>95%) and 
PHEV (>93%). 


FIGURE 1 Phylogenetic analysis by maximum likelihood method based on the partial RNA-dependent RNA polymerase (RDRP) gene (a), 
complete genomic sequence (b) and putative 4.8 kDa non-structural protein (c). (a) The phylogenetic tree was drawn based on partial RDRP 
nucleotide sequences (256 bases) of the three water deer coronaviruses, (indicated by black circles) and other reference sequences; (b) The 
phylogenetic tree was drawn based on genomic sequence of water deer coronavirus W17-18 and other reference sequences (sequence 
similarities enclosed in parentheses); (c) The phylogenetic tree was drawn based on amino acids sequences of the putative 4.8 kDa non- 
structural protein of water deer coronavirus W17-18 and other reference sequences 
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@ W17-18 CoV 
(a) HECV-4408 (FJ415324) 
821 @ W17-15 Cov 
BCoV(U00735) 
@ W17-12 CoV 
100 HCoV-OC43(NC 005147) 
38L— PHEV(NC 007732) 
op HCoV-HKU1(NC 006577) Betacoronavirus 
55 MHV(AB551247) 
o SARS CoV-Tor2 (JX163928) 
100 SARS-like Bat CoV-RA (00412042) 
57 33 Bat CoV-273 (DQ648856) 
BatCoV-HKUS(NC 009021) 
MERS-CoV(NC 019843) 
88 BatCoV-HKU4(NC 009019) 
100 r IBV(NC 001451) 
onavirus 
TCV(NC 010800) 
ear: TGEV(AJ271965) 
99 CCoV(KP981644) 
FCoV(EU186072) 
100 r BatCoV-1A(NC 010437) 
3 26 BatCoV-1B(NC 010436) 
BatCoV-HKU2(NC 009988) ; 
onavirus 
34 HCoV-NL63(DQ445911) 
85 23 HCoV-229E(JX503061) 
BatCoV-RoHKU10(NC 018871) 
` BatCoV-512(NC 009657) 
36 BatCoV-HKU6(0Q249224) 
23 PEDV(KC189944) 
—— 
0.1 
(b) 100) &F424621.1 Sable antelope coronavirus US/OH1/2003, complete genome (98.52%) 
4? EF424622.1 Giraffe coronavirus US/OH3-TC/2006, complete genome (98.48%) 
79 NC 003045.1 Bovine coronavirus, US/BCoV-ENT/1998, complete genome (98.56%) 
a EF424619.1 Bovine coronavirus US/E-AH187/2000, complete genome (98.79%) 
FJ425187.1 White-tailed deer coronavirus US/OH-WD470/1994, complete genome (98.81%) 
100 FJ938066.1 Bovine respiratory coronavirus bovine/US/OH-440-TC/1996, complete genome (98.69%) 
100 § F4425184.1 Waterbuck coronavirus US/OH-WD358-TC/1994, complete genome (98.67%) 
10! FJ425189.1 Sambar deer coronavirus US/OH-WD388/1994, complete genome (98.69%) 
ar Water deer coronavirus W17-18 
100 FJ415324.1 Human enteric coronavirus 4408, complete genome (98.19%) 
49t U00735.2 Bovine coronavirus US/Mebus/1972, complete genome (98.20%) 
100 KF530095.1 Human coronavirus OC43 strain OC43/human/USA/912-6/1991, complete genome (95.77%) 
100 L KX344031.1 Human coronavirus OC43 isolate LRTI 238, complete genome (95.52%) 
DQ011855.1 Porcine hemagglutinating encephalomyelitis virus strain VW572, complete genome (93.09%) 
NC 019843.3 Middle East respiratory syndrome coronavirus, complete genome (62.46%) 
100 NC 004718.3 SARS coronavirus, complete genome (61.83%) 
NC 002645.1 Human coronavirus 229E, complete genome (59.33%) 
weg. 
0.1 
(c) BCoV-R-AH187-ABP38288 






HECV-4408-ACJ35491 
BCoV-Bubalus/Italy/179-ABV 74056 
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FIGURE 2 Comparison of CoV genomes. The genome organization is schematically reported, grouping the viruses according to their genus 
or lineage. Grey-coloured boxes represent accessory open reading frames (ORFs) of the CoV genomes 


The full genome sequence-based phylogenetic analysis by maxi- 
mum likelihood method revealed that the virus was closely related 
to BCoV, bovine-like coronavirus and HECV-4408 (Figure 1b). 
Distinctively, the N2 protein (208 amino acids), which was solely 
found in the BCoVs, was detected in the sequenced genome, 
whereas the predicted 4.8-kDa non-structural protein in water deer 
coronavirus strain W17-18 was elongated due to a 13-nucleotide 
deletion compared to other BCoVs (Figure 2). Therefore, the phylo- 
genetic analysis based on the 4.8 kDa- or 4.8 kDa-like non-structural 
protein amino acid sequences of other reference viruses infers that 
water deer coronavirus W17-18 has created a unique clade within 
the BCoVs (Figure 1c). Collectively, the newly found water deer 
coronavirus in this study was closely related to BCoV showing minor 
unique genetic characteristic in the 4.8 kDa-like non-structural 
protein. 

As the water deer coronavirus was detected from nasal swabs 
during the winter season, it could be assumed that the water deer 
coronaviruses might also cause respiratory disease and circulate dur- 
ing the winter season like other BCoVs. The majority of BCoV cases 
previously detected in wild ruminants have been in captive animals. 
In this study, BCoV was detected in non-captive wild water deer, 
which may provide further evidence that BCoV could be transmitted 
between freely ranging wild ruminants and farm ruminants. These 
findings may attract public attention for BCoV transmission between 
wild water deer and farms during winter season. In addition, BCoV- 
like human coronavirus (HECV-4408) was demonstrated to be able 
to infect gnotobiotic calves and provide cross-protection against 
BCoV infection (Han, Cheon, Zhang, & Saif, 2006), and interspecies 


transmission of BCoV in turkey poults was also demonstrated under 


experimental condition (Ismail, Cho, Ward, Saif, & Saif, 2001). There- 
fore, interspecies transmission ecology of BCoV among livestock, 
wild animals and humans should be further investigated to control 
the virus. 

We conclude that the water deer coronavirus is a new strain of 
BCoV, rather than a new species. This is the first report of BCoV 
detection in non-captive wild water deer and presents a possible res- 
piratory infection risk to water deer in Korea. Water deer (H. inermis) 
have recently been recognized as the most common wildlife to 
approach farmhouses and livestock barns in Korea (Kim et al., 2014). 
BCoV transmission between water deer and livestock farms should 
not be overlooked. Therefore, effective biosecurity systems in 
bovine farms should be investigated to prevent interspecies trans- 


mission of BCoV between farm ruminants and wild water deer. 
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